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Abstract 

We carry out simulations of the collision of two components of an adiabatically divided, quasi-2D 
EEC. We identify under, over and critically damped regimes in the dipole oscillations of the com- 
ponents according to the balance of internal and centre-of-mass (cm.) energies of the components 
and investigate the creation of internal excitations. We distinguish the behaviour of this system 
from previous studies of quasi- ID EEC's. In particular we note that the nature of the internal 
excitations is only essentially sensitive to an initial phase difference between the components in the 
over damped regime. 
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I. INTRODUCTION 



Quantum interference (QI) in atomic BECs has been an important topic since the exper- 
iment carried out in Ref.[ ] demonstrated the matter- wave nature of condensates. Interfer- 
ence effects have been used to investigate decoherence [2] and vortex production [3, 4] while 
underpinning many potential applications in metrology, gravitational, temporal and rota- 
tional sensing [5-7]. Most experiments on atomic interference have used freely expanding 
atom clouds [8-10]. An alternative approach, however, involves controlling the the motion 
of atom clouds using time dependent trapping potentials [11], often generated by atom chips 
[12], which can be used to study atom-surface interactions and Josephson phenomena[13]. 
Many discussions of confined-atom interference consider the noninteracting situation [14, 15] 
or do not address the deleterious effects of interatomic interactions, including dephasing, 
collisional shifts, and phase diffusion [16-22]. However, in Ref.[ ] it was shown that the 
recombination process of a split condensate is very sensitive to atomic interactions. In this 
paper we will study the damping of dipole oscillations of colliding, coherent, condensate 
clouds and the consequent creation of internal excitations in the clouds. 

Quasi-2D EEC's can be formed in a highly anisotropic trap of the form 

Krap (r)=y^(x,2/)+K(z) 

= ^mul {x^ + y^) + ^mulz^ , (1) 

with (jUz ^ In this limit the dependence on the z-coordinate is quenched such that 
the order parameter has the form ^!{x,y,z]t) = a/]V0(^, ?/; exp{— z^/2/^}/(7r/^)-^/^ 

where (f) is normalized to unity. The dynamics of such a 2D BEC is governed by the CP 
equation 



where g2 = g-iD/V^h = 2\/2ti a si z^jJz is the 2D coupling constant, 1^ = y^h/mujz is the 
oscillator length in the tightly confined direction and is the bulk 5-wave scattering length. 

If a BEC is prepared in such an anisotropic trap by cooling to its ground state and the 
trap is then suddenly displaced through a distance —A in the x-direction then the cloud 
will execute simple harmonic motion without any internal distortion. In the absence of 
interactions the initial state is a simple gaussian which evolves according to 

0(x,|/;t) = 0o(a;-X(t),?/)e^^«^/'^ (3) 

where 

M^,y) = \A2^M-{^^ + y^)/^^l} , (4) 



where l\ = h/mu± and 



X{t) = Acos{iJi_t) (5) 
P{t) = mu±A sin {co±t) (6) 



which are the position and momentum of a classical particle of mass m prepared at rest at 
position a; = A in the same potential. The inclusion of interactions only alters the form of 
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the initial state, (po^x^y) but does not alter the forms for X{t) and P{t) and, hence, does 
not qualitatively change the nature of the subsequent dipole oscillation (see section 12.2 of 
[24]). 



We are considering the situation in which a quasi-2D BEC, previously prepared in an 
anisotropic trap, has been adiabatically divided into two spatially disjoint components while 
retaining its overall coherence. The two components are supposed to experience slightly 
different chemical potentials so that they acquire a phase difference proportional to the 
time, th^ for which they are held apart: 9 = Sfit. The trap is then suddenly switched 
back to a simple parabolic trap. This is equivalent to preparing a BEC in a state that is a 
superposition of two single clouds (as above) displaced by equal and opposite distances ±A 
and with a phase difference, 9 

$(x, y;t = 0) = ^ ((/)o(x + A, ?/) + e^Vo(x - A, y)) . (7) 

The system then evolves according to the GP equation. The case of a non-interacting 
condensate is, of course, exactly solvable as the two components evolve independently. At 
time t > the state of the non-interacting system can be written as 

y,t) = ^ (e-^^«^/Vo(a; + X{t),y) + e^^e^^«^/Vo(a; - X{t), y)) (8) 
with X{t) and P{t) as given above. The corresponding density is then 

n(x, y, z] t) = 712 {x, y] t) e^-^'/^'VV^ (9) 

where 



n2{x,y]t) = A^|$ 
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= ^e-(^"«'«'<")/'i {cosh [^) + cos (H^ 4) } . (10) 

The two components are maximally overlapped at times = + |) ttcj^ when 

nfe,;T„) = A'^e-('^+»')"i|l+cos(?^ + 0} . (11) 

Corresponding to perfect rectilinear fringes in the x-direction with spacing given by the de 
Broglie relation L = h/2m(jj^/\. The initial relative phase difference simply shifts the fringe 
pattern. As is well known the shot-to-shot reproducibility of the fringe positions is a test of 
the coherence of the two components. 



If we consider the energetics of the non-interacting case we can express the mean-field 
energy of the initial state as Et = Ek + Ep where 

= ^ / ^'"1^ l^'^l' = ^2^ + ^\^^ ^^^^ 
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is the kinetic energy and 

Ep = N I dh-mujl |r|^ |0|^ = N-mujl {X{t)f + A^-^^ (13) 
J 2 2 2 

the (trap) potential energy. We can regroup these into two conserved quantities 

Ecm = N^mcvlA^ (14) 
is the energy associated with the bulk motion of the two clouds and 

^int = Nhw^ = (15) 

is the internal energy of the two clouds (in this case purely the confinement energy). 

The inclusion of interactions between the atoms will, naturally modify this picture. While 
a single BEC, when subject to a displaced trap potential, will execute undamped simple- 
harmonic motion with no change of shape (dipole mode) a pair of separated components 
will only only undergo independent dipole oscillations until they spatially overlap. The 
subsequent dynamics will be more or less complex depending on the relative magnitudes of 
the centre of mass (cm.) energy and the mutual interactions. 

For moderate interaction strengths we would expect the initial components to be well- 
described by the Thomas-Fermi approximation 

^ [R^, - |r| ) r < Rtf (^g) 
|r| > Rtf 

where, as always, Rtf is chosen to give the correct total number of particles which yields 

Rtf = V2/^7 (17) 
where the 2D Thomas-Fermi parameter is 

7^(^V' . (18) 



Hence the conserved energy will be Et = Ecm + ^int where 

^cm = N^mujlA^ (19) 
^int = Etf = ^f^TFN = Nhw^^^ (20) 

These will not be separately conserved and we expect an irreversible transfer of energy 
from Ecm to E^^t leading to damping of the dipole oscillations of the two components at the 
expense of the creation of internal excitations within each cloud. 

We seek to define an effective measure of the relative importance of interactions which 
can be directly evaluated in a simulation. We define a parameter, 77^, at each turning point 
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Tn = rni/uj^ of the dipole oscillation (at which the cm. kinetic energy is zero) as 

Ep lmu;i/|rp|(/>(r,r,)|'D2r ' ^ ' 

In the non-interacting limit {§2 = 0) this is identically zero. In the Thomas Fermi regime it 
will be 

— 1 — 22 

^muiAl 

where is the displacement of the cm. of the cloud from the centre of the trap at time 
Tn and 6ex is the excitation energy per atom of the cloud (i.e. the excess internal energy of 
the cloud over the TF energy per particle) . Since the scattering will act to damp the dipole 
oscillations (decreasing A) and create excitations within the cloud (increasing cex) Vn will 
increase with n. 

For low density/ weak interactions {r]o <^ 1) one expects the two components to undergo 
underdamped oscillations and to exhibit behaviour which is qualitatively similar to the 
non-interacting case. In particular we would expect curved fringes while the two clouds 
are passing through one-another. The scattering processes arising from interactions will 
produce a frictional force on the cloud tending to reduce its velocity and hence increase 
the fringe spacing. The frictional force will be largest in the centre of the trap (close to 
y = 0) where the density is highest. In addition, the interaction will lead to an over-pressure 
in the fringe maxima which will tend to fatten the fringes in the trap centre. At higher 
densities/stronger interactions we expect to find an overdamped regime in which the two 
clouds collide and merge into a single cloud containing a variety of excitations seeded by the 
fringe pattern associated with the opposite pre-collision velocities of the clouds. As we will 
see it is in this overdamped regime that we observe significant qualitative sensitivity to the 
initial phase diflFerence, 9^ between the two components. A recent experiment [25] showed 
that the recombination of a split condensate held by an atom chip leads to a heating of the 
atomic cloud which depends on the relative phase of the two components. This qualitatively 
agrees with the analytical prediction that the recombination process can lead to exponential 
growth of unstable modes [23]. This efltect, however, is not obvious in the underdamped 
regime. The role of the relative phase can make the dynamics of the interferometer very 
different in the overdamped regime [25, 26]. 

In this paper, we present detailed numerical calculations, which now enable us accurately 
to characterise the behaviour of the interferometer as a function of the distance and the 
relative phase between the two disk-shaped condensate components. We carry out the sim- 
ulations in the underdamped, overdamped and critically damped regimes. Our key findings 
are as follows. 

Underdamped case The cm kinetic energy is much larger than the interatomic interaction 
energy. In this case the system develops distorted rectilinear fringes. 

Critically damped case The cm kinetic energy and the interatomic interaction energy are of 
the same order. This means that nonlinearity is important and the effect of scattering 
becomes noticeable making the dynamics of the system complicated. After the clouds 
interfere for the first time, topological excitations are created in the system. In both 
the underdamped and critically damped cases we find that the initial relative phase 
does not affect the system qualitatively. 
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FIG. 1. The equipotential lines of a 2D double well potential with narrow cental barrier along 
about X = 0. The parameters used are = 0.1, Oyjl^ — 10, U^^jhuj^ — 50. 



: Overdamped case The cm energy is much smaller than the interaction energy, here the 
system undergoes a merging process and is much more sensitive to the initial relative 
phase. For zero initial relative phase the two components can merge into a single cloud 
containing collective excitations, while for a tt initial relative phase the merging leads 
to topological excitations. 



II. SYSTEM AND METHODOLOGY 



In experiments, a trapped condensate can be split into two components adiabatically, for 
example by using a tailored magnetic potential from an atom chip [27], by the introduction 
of an optical barrier via a shaped blue-detuned laser [ ] or by passing counterpropagating 
red-detuned laser beams through an acousto-optic modulator driven at RF frequency [25] 
to form a double well potential in the horizontal plane. The two components of the split 
condensate are held in an asymmetric double well potential for varying hold times to set the 
relative phase between them. The recombination of the two components after splitting can 
be realised by reducing the RF frequency for variable recombination times. 

We based our simulations on a condensate of = 2 x 10^ Rh atoms in a trap potential 
with cj^ = 27r X 5i7^ and cj^ = 27r x lOOi/^. To simulate the underdamped and critically 
damped cases, we obtained the configuration, (\)^[x^y\ of the order parameter which min- 
imised the mean field energy subject to the constraint J |(/)o|^ = 1. The initial state for 
the simulations was then 

^; t = 0) = -L (0o(x + A, ?/) + e^Vo(^ - A, y)) . (23) 

The order parameter was then allowed to evolve in the quasi-2D trap potential. To reach 
the overdamped regime required A < i?TF so another method of preparing the initial state 
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FIG. 2. Normalised density profiles of the two condensate components without initial phase differ- 
ence colliding at the bottom of the trap, (a) noninteracting case with A = 7.5/^; (b) interacting 
case with A = 7.5/^; (c) interacting case with A = 4/^. 

was used. The lowest energy configuration of the order parameter was found in a potential 

1 _x2_^2 

VDw{x,y) = -mujl{x^ + y^) + U^e ^ ^. (24) 

where the second term represents the effect of a narrowly focussed, blue-detuned laser applied 
to split the condensate. The height of the central barrier is given by the amplitude of the 
Gaussian function [/q, and the distance between the two wells can be adjusted using the 
parameters and Gy. We can see two distinct potential minima clearly in Fig. 1, which 
have the same shape as the RF-induced potential in the experiment in Ref.[2^j]. The phase 
difference between the two parts can be achieved by using phase imprinting method [30] in 
experiments. The parameters we used in our simulations are (Jx/l^ = 0.1, cfy/l^ = 10, and 
Uo/hu)j_ = 50. We then used this configuration as the initial state for real-time evolution of 
the GP equation, representing the sudden removal of the laser. 

In our discussion we will always keep the interaction strength and the number of atoms 
in the system fixed. The only quantities we change are the initial displacement and the 
initial relative phase of the condensate components. The ground state wave function of the 
condensate is obtained by the imaginary-time-evolution method [31] for the time-dependent 
GP equation (2). We use a reduced scattering length for ^^Rb^ = 0.3 x 5.4nm to see the 
structure of vortex excitation clearly in our simulation. To consider the effect of the relative 
phase on the dynamics, we keep the initial phase of one component fixed and adjust the 
initial phase of the other from to 27r. 



III. NUMERICAL SOLUTION AND RESULTS 

In the numerical investigations we used l± as the unit of length and 1/(jJ± as the unit of 
time. Hence we define dimensionless wavevectors k = kl±^ times i = (jj±t and frequencies 

UJ (jU/(jUj_. 



7 





FIG. 3. Normalised density profiles of the colliding components in the trap without initial phase 
difference. The initial separation A is 7.5/^ and 4/^ (b)(d) respectively. Times are (a) 

i = 2.15, (b) i = 3.05 and (c)f = 4.25, (d) t = 3.5. 



A. Scattering effect in the under damped regime 

As we have discussed, the importance of the scattering between atoms depends on the 
competition between the cloud kinetic energy and the interaction energy, which in turn affect 
the dynamics of the system. 

In Fig. 2, we show the density profiles at the point of maximum overlap {t = 7r/2) of 
the two components with zero initial phase difference starting from different displacement 
positions. The bright region is a low density area, while the dark region is a high density 
area. For a noninteracting system there is only kinetic energy so 77 = 0. In this linear regime 
the unperturbed condensates always show straight fringes during the interference as seen 
in Fig. 2(a). For interacting systems when the displacement is large, the clouds will have 
maximum cm. kinetic energy, i.e. high relative peak velocity when they begin to interfere. 
The system still responds linearly. The duration of the collision is correspondingly short and 
the effect of scattering is very small. In Fig. 2(b), the initial separation is set to be A = 7.5/^. 
We can see that interference fringes are nearly straight. With smaller initial separations, 
the components reach a smaller peak velocity, the interaction energy is comparable to the 
kinetic energy and the duration of the collision is correspondingly longer. As is shown in 
Fig. 2(c), where A = 4/^, the distortion of the fringes has become quite sizable. The 
width of the fringes increases with the nonlinearity, while the peak density does not change 
dramatically in comparison to that of the fringes in the linear situation (Fig. 2(b)). The 
nonlinear situation however gives us fewer fringes associated with the damping of the dipole 
oscillation of each cloud. The expansion rate of the central fringe is highly nonuniform along 
^-axis due to the non-uniform density distribution. At the centre of the clouds, where the 
density is highest, the fringes expand most. The spacing of fringes are close to their non- 
interacting value towards the edge of the cloud (low density region) where the interaction 
is negligible. This leads to the central fringe taking a lenticular shape, accompanied by the 
bending of the other fringes. There is only the central fringe left eventually before the two 
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FIG. 4. The normalised density profiles (the left panels) and the corresponding phase plots (the 
right panels) of the colliding condensate components in a trap during the second recombination. 
This is the continuation of the simulation shown in Fig. 3. From (a) and (b) we can see clearly 
the formation of two vortex pairs located at the left and the right side of the centre of the trap, 
curved solitons arise in the centre of the trap. These vortex pairs approach each other until the 
vortices with opposite charges annihilate, forming a soliton ring as seen in (c) and (d). The vortex 
ring will decay into vortices thereafter. Times are (a,b) t — 4.15, (c,d) t — 4.5. 

components start to separate. 

In Fig. 3 we show the density profiles of the two components without initial phase 
difference after their first collision, where the initial separation A is 7.5/^ and 4/^ 

(b)(d) respectively. The effect of interactions is not obvious in the linear regime (A = 7.5/^). 
As seen in Fig. 3(a) the condensate expands a little bit in the ^-direction which means that 
some atoms gain a y-component to their momentum via scattering. This momentum is 
small, however, due to the short collision time, which means that the scattering induced 
distortion of condensates is not obvious after the first collision. After the first collision the 
clouds pass through each other; each moving towards the original position of the other. The 
density profile of the condensate will vary in a manner similar to the the breathing mode 
and the components will undergo a damped cm. oscillation. Fig. 3(c) shows that the clouds 
regain their original shape with a small distortion at t = 4.25 before their second collision. 

When the initial separation is reduced, the maximum cm. kinetic energy is correspond- 
ingly smaller and the nonlinear interaction term in Eq.(2) becomes more important. We 
reach the critically damped regime (e.g., A = 4/^) where the cloud energy is comparable to 
the interaction energy. As seen in Fig. 3(b) the modulation of the density in the middle of 
the component (high density area) moves rapidly in the ^-direction due to strong scattering, 
leading to two density peaks appearing symmetrically about the x-axis for each component 
(see Fig. 3(b)). Due to the trap potential, self-interference occurs when these two density 
maxima move back towards one other as seen in Fig. 3(d). The most distinct feature in 
this regime is that the clouds no longer pass cleanly through each other and regain their 
original shape after their collision. We can see clearly that there are more structures in 
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FIG. 5. Oscillations of condensate atom number Ny for different initial conditions. 

the condensate in Fig. 3(d) than in Fig. 3(c). The interference fringes will be completely 
distorted when the second collision begins and this leads to excitations in the form of vortex 
pairs and soliton rings as shown in Fig. 4. 

In Fig. 5, we show quantitatively the scattering effect by plotting the dynamics of Ny^ 
which is the number of atoms with ^/-component of momentum \ky\ > lOdk where 27r/A; is 
the spatial resolution of the simulation. As seen in the figure, the initial value of ky is 
because the condensate components are set to move along the x-axis. The black solid line 
is the absolute value of a sinusoidal curve fitted with fixed frequency uj = 1 which indicated 
that there are two collisions during one dipole period. For all three distinct initial conditions 
there is only one peak during the first half period because the scattering occurs after the 
two components collide at the bottom of the trap (t ^ 7r/2). The peak value increases with 
decreasing A. After the first half period two peaks appear. The first one is induced by 
the self-interference in the y-direction of each component, and the second one arises from 
the collision of the two components. For the underdamped case (the line marked with open 
circles) the self-interference induced scattering is always weaker than the collision induced 
scattering before the conditions of the linear regime are violated (recall that rjn increases 
with n), while it is true for the critically damped case (the line marked with cross and 
the line marked with dot) only in the second half period. In this intermediate regime the 
creation of solitons and vortices after the second half period suppress the scattering along 
the x-direction (collision induced scattering), which makes the second peak lower than the 
first one. 

Generally in the above regimes the initial phase difference of two components will not 
change the physics intrinsically. Especially for the underdamped situation, only the location 
of the fringes varies with the relative phase while the period and spacing of the fringes 
are not changed. In the underdamped regime, the trajectory is independent of A9 indicating 
that the initial relative phase will not change the scattering properties at all. In the critically 
damped regime the only difference between the = (dotted line) and A9 = n trajectories 
is the peak values. An interesting point is that the peak values of A9 = tt is always smaller 
than those of A9 = 0, which is different from the results for the merging process (overdamped 
collision) in Ref.['^ ] where the atom loss after recombination increased as the relative phase 
increased from to tt. 

If we reduce the distance between two components further, the interaction energy will 
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FIG. 6. The normalised density profiles and phase plots of the condensates prepared in a double 
well potential which are released into a harmonic potential well at time t = 0. The initial relative 
phase is A6> = 0. The left panels are the density profiles of the condensates prepared in a double 
well potential (a) and then released in the harmonic potential (c, e, g, i). (d, f, h, j) are the 
corresponding phase plots of (c, e, g, i) respectively, (b) is the density profile of the ground state 
condensate with 2N atoms in the harmonic potential well. The density profile of the dynamic state 
(g) is close to the ground state (b), but with some distortion and a few excited atoms around it. 
Density and phase modulations are obvious in plots (e) and (j). Times are (a) t = 0, (c, d) t — 0.5, 
(e,f) 1.0, (g, h) t = 1.57, (i, j) t = 10. 



dominate the system. In this overdamped situation the dynamics of the system could be 
very different from the cases we have considered so far. The two components could merge 
into one condensate after the first collision. In the next section we will look at this nonlinear 
regime where the initial relative phase plays a more important role. 



B. Phase-sensitivity of excitations in overdamped collisions 

In Ref.[25, 26], it was reported that for elongated components that are tightly confined in 
the radial directions, when the central barrier potential is suddenly turned off, two compo- 
nents merge into a single cloud. When the initial relative phase between the two components 
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is zero, the merged cloud is in its ground state with respect to the trapping potential. When 
the initial relative phase is tt, the merged condensate contains a dark soliton. 

For 2D disk-shaped clouds, the size of the condensate in the merging direction (x-axis 
in our case) is large. As we have seen, when the separation between the cm. of the two 
components is comparable to the diameter of one cloud, i.e. A ~ Rtf^ the cloud kinetic 
energy that the components will have is comparable to the interaction energy. This energy 
is high enough to prevent the two components merging into one condensate after the first 
collision, as they will not lose sufficient momentum in the merging direction due to scattering. 
If we want to investigate the merging of 2D condensate components we need to enter the 
overdamped regime where the interaction energy dominates the system, i.e. A < Rtf- This 
means that there will be an overlap between the two components which depends on the 
height and width of the central barrier of the trap. Vow (see equation 24). 

As we have seen, solitons are created during the interference process. We know that 
soliton stability depends on the nonlinearity and geometry of the medium [ ]. In general, 
solitons in BECs are always thermodynamically unstable, and dark solitons in BECs are 
also expected to be inherently dynamically unstable [34-37]. The velocity of a soliton Vg is 
[38] 

Vs/vc = cos{A9/2) = [1 - Ud/n]'/^. (25) 

'^c = (5'2D^/^)'^^^ is the speed of sound, n = |(/)op is the unperturbed condensate density, 
and rid is the difference between n and the minimum density in the centre of the soliton. For 
a given density, increasing the relative phase from to tt will decrease the velocity of the 
dark soliton. The two limiting cases are ^ = and = tt, where the velocity of a soliton is Vc 
and respectively. The formation of vortex pairs inside the combined condensate depends 
strongly on the velocity of the soliton Vs and the density distribution of condensates. 
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FIG. 8. Sequence showing density profiles (left panels) and the corresponding phase diagrams 
(right panels) for components prepared in a double well potential undergoing a sudden merger in 
a harmonic trap. The initial relative phase is A6> = tt. Times are (a,b) I = 0.5, (c,d) I = 7.5, (e,f) 
f=8, (g,h) t = 15.2. 

Firstly, we consider the situation for which the initial relative phase between the com- 
ponents is 0. As shown in Fig. 6(a), there is an overlap between the components in the 
left well and that in the right well. The dynamics of the two components will lead to them 
interfering and merging after the central barrier is turned off. The spacing of the central 
fringe increases with time until the two components reach their maximum overlap and the 
shape of the central fringe is close to that of the ground state of a condensate with 2N 
atoms in the harmonic trap (see Fig. 6(g) and (b)). The outer fringes nearly vanish and 
the two components merge into one with some excited atoms being found near the edge of 
the cloud (see Fig. 6(g)). The scattering is sufficient due to the long collision time which 
ensures that the combined condensate can not separate any more. The expansion of the 
condensate is not as obvious as in Fig. 3(b) due to the small maximum cm. kinetic en- 
ergy of the components and the fact that the interatomic interaction dominates the system. 
However the cloud cannot reach its ground state, as would happen in the ID case, because 
the trapping frequency in the merging direction is too low. As predicted by Eq.25 the dark 
solitons move at and decay rapidly so the system should be dominated by density wave 
excitations. In Fig. 6(g) (h) we see that a state with density and phase modulations, but 
without topological defects, is obtained, meaning the system is in a collective excited mode. 

Secondly, we consider the situation where the initial relative phase is tt, but with the 
other conditions are kept unchanged. We find this merging process is completely different 
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FIG. 9. Oscillations of condensate atom number (red dotted line), Ny (green dashed line) and 
Nxy (thick solid line) for different initial relative phases. The cyan thin lines are sinusoidal curves 
fitted with frequencies a; = 1.66 (a), 1.6 (b), 1.3 (c), 1.0 (d), 0.6 (e). 

to the case where = 0, and that it is also distinct from the ID merging of condensates 
with relative phase tt reported by Ref.[ ] and Ref.[26], where the merging time is too fast 
(merging occurs in the tight confining direction) to make the soliton decay. In 2D cases the 
radial frequency of the trap is weak, and so will not quench mechanical excitations in either 
X- or ^-direction. As seen in Fig. 7(a) and (b), the condensates with relative phase tt are 
initially prepared in a double well potential. At time t = 0, the central gaussian barrier 
is turned off suddenly. For interference with tt relative phase there should be two central 
fringes. Due to the small initial separation of the two components, these two central fringes 
are very fat and nearly keep the original shape of the initial clouds. Our simulations show 
that two density waves propagate in opposite directions at a reduced amplitude as seen in 
Fig. 7(c). The initial soliton will eventually reach a point where its central density is zero (see 
Fig. 7(c) and (d)), and the phase difference is tt between the two parts, where a completely 
dark soliton (black soliton) is created. The velocity of the black soliton is zero. In Fig. 8, 
we give the snapshots of the time evolution of the condensates which show that the black 
soliton is not as stable as in the merging of elongated components. After approximately 
143ms (t 4.5) of real-time propagation, the black soliton begins to undulate. It is bent 
at first, then divides into two from the middle with atoms passing through the centre and 
propagating from left to right. The atoms in the right cloud are pushed by these additional 
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atoms and then two paths from right to left are formed, one above and one below the centre 
of the cloud, leading to the decay of the solitons via the snake instability [39, 40]. A pair 
of vortices form near the centre of the condensate and subsequently the other two vortices 
nucleate at the edge of the cloud, move around the edge to the left and towards each other, 
and then annihilate at y = 0. We note that due to the initial tt relative phase the two 
components are still distinguishable by the phase difference as shown in Fig. 8(f) and (h). 

For initial relative phases between and tt, the initial dark soliton propagates in the 
positive X direction, and curves with the expansion of the central interference fringe. This 
curvature is similar to that observed in Ref.[38]. However the dark soliton will eventually 
decay into vortex pairs instead of oscillating in the trap. 

The nonuniform density distribution in our system makes the speed of a soliton a function 
of which leads to the curvature of the solitons. the formation of vortex pairs inside the 
combined condensate would require longer lifetimes (i.e. smaller velocity) for the solitons. 

In Fig. 9 Nx is the number of atoms with > lOd/c, and Nxy is the number of atoms with 



and 



larger than lOdk. We can see that Ny is almost unchanged until the soliton 

3egan to decay (green dashed line), which means that the scattering in the y— direction is 
suppressed by the formation of solitons. The lifetime of solitons increases as is increased 
from to 271 as predicted by Eq.25. When is smaller than 7r/2 the oscillation frequencies 
of Nx and Ny are nearly the same after the decay of solitons as is that of the A^^^^ oscillation. 
When A9 = the velocities of the accompanying density waves are close to Vc [41] and 
the corresponding frequencies are close to 2uj±. As shown in Fig. 9(d) and (e), when A^ is 
above 7r/2 the formation of vortices inside the combined condensate as seen in Fig. 8 (g) will 
diminish the amplitude of the number oscillation, which means that the presence of vortex 
pairs in the combined condensate suppresses density waves. 



IV. CONCLUSION 

In this paper we have discussed the excitations created when a quasi 2D BEC is first 
divided into two coherent components and which are then driven to collide by a harmonic 
potential. The collision and merging processes are quite different to those in ID or quasi- 
ID (elongated) systems. However the fact that the dynamical excitations depend not only 
on the initial separation of the components but also on their initial relative phase is the 
same as the ID case reported in Ref.[25, 26]. Generally the dynamics of the system are 
decided by the competition between the cm. energy and the internal energy of the colliding 
components, but it is not simply that the smaller is the maximal cm. kinetic energy, the 
easier it is to excite vortices in a system of interfering BECs. When the kinetic energy dom- 
inates, the scattering is less effective due to the short collision time and linear interference 
is observed. There are no dynamical excitations after the two components are separated 
and the excitation properties are not phase sensitive. If the maximal cm. kinetic energy is 
comparable to the interaction energy, the scattering is effective which makes the transverse 
self-interference of each component obvious. The combination of the interferences in both 
X— and direction makes the dynamical excitations in the system quite complicated. How- 
ever the role of initial phase difference is still limited to the position of fringes. When the 
interaction energy is much larger than the kinetic energy (merging process) the scattering 
is still less effective because of the small maximum cm. kinetic energy of the components 
and the longer lifetime of solitons. Nxy is smaller than 4% during our simulations. The 
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two components can merge into one in this case. The final states are, however, excited 
in diflFerent ways, depending on the initial relative phase. We found that, only when the 
initial relative phase is larger than 7r/2, will there be a vortex pair which can persist in the 
combined condensate. Otherwise, all the vortex pairs formed in the merging process will 
only move around the edge of the combined condensate and escape so that the system is left 
with density wave excitations. The formation of vortices inside the condensate suppress the 
density wave oscillations. 
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